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We study the zero temperature ground state of a two-dimensional atomic Fermi gas with chemical 
potential and population imbalance in the mean-field approximation. All calculations are performed 
in terms of the two-body binding energy e_g , whose variation allows to investigate the evolution from 
the BEC to the BCS regimes. By means of analytical and exact expressions we show that, similarly 
to what is found in three dimensions, at fixed chemical potentials, BCS is the ground state until 
the critical imbalance h c after which there is a first-order phase transition to the normal state. 
We also show that for a fixed ratio eb/ef, where cf is the two-dimensional Fermi energy, as the 
density imbalance m is increased from zero, the ground state evolves from BCS to phase separation 
to the normal state. At the critical imbalance m c phase separation is not supported and the normal 
phase is energetically preferable. The BCS-BEC crossover is discussed in balanced and imbalanced 
configurations. Possible pictures of what may be found experimentally in these systems are also 
shown. 
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I. INTRODUCTION 

In the last few years, great experimental advances have 
permitted the manipulation of trapped neutral ultracold 
two-spin-components atomic Fermi gases with tremen- 
dous accuracy. The essential technics under full domain 
are the cooling, trapping, the control of the number 
of atoms in the sample, and the tunning of the inter- 
atomic (s-wave) interactions via the application of an ex- 
ternal uniform magnetic (Feshbach resonance) field 
This allowed the investigation of the crossover from the 
Bardeen-Cooper-Schrieffer (BCS) phase of weakly bound 
Cooper pairs to the Bose-Einstein condensate (BEC) 
phase of strongly bound diatomic molecules in three- 
dimensional (3D) trapped Fermi gases 

These many-body quantum gases are even more ex- 
citing when the two-components have mass, chemi- 
cal potential or population imbalance, since this situ- 
ation can be met in many areas of physics, from con- 
densed matter to high-density quark matter. Many ex- 
otic phases have been proposed as the ground state of 
these imbalanced Fermi gases, such as the homogeneous 
Sarma Q or breached-pair state @, the Fulde-Ferrel- 
Larkin-Ovchinnikov (FFLO) state with modulated order 
parameter Q, phase separation of superfluid and nor- 
mal c omp onents in real space Q , deformed Fermi sur- 
faces and mag netized superfluid in the BEC side of 
the resonance [ll[. Recent experiments in imbalanced 
atomic Fermi gases confined in 3D harmonic traps have 
observed that these systems phase separate into a un- 
polarized superfluid core surrounded by a normal outer 
region [12h16| |. 

Another very important parameter under control is the 
trap geometry, which permits the experimental investi- 
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gation of fcrmionic atoms in lower dimensions. A two- 
dimensional (2D) trapped Fermi gas may be achieved 
in experiments by flattening magnetic or dipolar con- 
finements [l7| . by trapping atoms in specially designed 
pancake potentials jl8j| , b y radio frequency (rf)-induced 
two-dimensional traps [19(, or by gravito-optical surface 
traps [13]. 

Motivated by these experimental progresses, some re- 
cent pa per s have investigated two-component Fermi gases 
in 2D |2l| - [24| . One-dimentional (ID) two-component 
Fermi systems have also received attention recently both 
theoretically [25|, [26[ and experimentally [27j . Besides, 
spin imbalance has also been investigated in the con- 
text of metallic systems as, for example, electron-hole 
bilayer [H, [2^ and multi-band [3(| systems. 

In this paper we study the zero temperature (T) 
ground state of a 2D atomic Fermi gas with chemical 
potential and population imbalance in the mean-field 
approximation. We show that in the grand canonical 
ensemble the fundamental state is BCS up to a criti- 
cal value of the chemical potentials asymmetry h c , after 
which there is a quantum phase transition to the nor- 
mal state. In the canonical ensemble we study in detail 
the issue of phase separation. We show that the nor- 
mal phase is always unstable to phase separation when 
the density asymmetry m is less than the critical im- 
balance m c . All calculations are performed in terms of 
the two-body binding energy e^, whose variation allows 
to investigate the evolution from the BEC to the BCS 
regimes. Since the recent 3D trap experiments of Shin 
ct al. [I?! have not observed the FFLO state, we do not 
consider in this work the possibility of this phase in 2D. 
The observation of the FFLO state is indeed an experi- 
mental challenging, since it is predicted to exist only in 
a narrow window of asymmetry between the particle's 
chemical potentials. Nevertheless it is worth to mention 
that significant FFLO-like correlations are expected to be 
present in ID systems for a range of polarizations (27j . 
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The main results of this paper may be summarized as: 

I. We show that in a zero T unbalanced 2D gas of 
fcrmionic atoms, BCS is the ground state up to chemi- 
cal potentials imbalance h c at which there is a first-order 
phase transition to the normal state. 

II. We demonstrate that, as found theoretically and ex- 
perimentally in 3D, for a fixed two-body binding energy 
e#, as the imbalance m is increased from zero the stable 
states are BCS, phase separation and normal phase. 

III. For balanced systems, the BCS-BEC crossover is 
governed only by e^- However, for unbalanced systems, 
the BCS-BEC crossover will depend on es and the im- 
balance m. For a fix value of es, the crossover is possi- 
ble only for asymmetries supported by phase separation. 
Above a critical value of the number densities imbalance 
m c , phase separation is destroyed (in favor of the normal 
phase) as well as the path to BEC regime. 

The paper is organized as follows: In Sec. [II] we present 
the model Hamiltonian and derive an analytical expres- 
sion for the zero temperature thermodynamical potential 
in the mean-field approximation for the two-component 
gas of fermionic atoms in 2D. In Sec. Mil wc employ 
this expression to investigate the stability of the possi- 
ble phases of 2D imbalanced Fermi gases, both at fixed 
chemical potentials and densities. We also discuss the 
BCS-BEC crossover in balanced and imbalanced systems. 
Sec.lTVlis devoted to the conclusions. 



II. THE MODEL HAMILTONIAN AND THE 
MEAN-FIELD THEORY 

We start by considering a 2D nonrelativistic dilute sys- 
tem of fermionic atoms of mass M , with two hypcrfinc 
states labeled as a =t,|- This spin \ and 1 mixture 
could, in principle, be done with the two lowest hyperfine 
states of 6 Li atoms, as in the 3D experiments fl3l. Ifil. [l6j . 
Their single-particle dispersion relations are given by 

h 2 k 2 

£ fc = T/TT' Throughout the paper we set h = 1. Since 
at low densities (appropriate to describe the ultracold 
trapped Fermi gases we are interested) the form of the 
potential is not probed, it can be considered that the 
atoms interact via a contact interaction which is mod- 
eled by the following pairing Hamiltonian: 



as n a = fi + ah, where ah = ±h. Then, the chemical 
potential \i a fix the number densities n CT of the differ- 
ent fermions. The new dispersions for the free species a, 
relative to their Fermi energies, are e£ = £fc — The 
interaction Hamiltonian is given by 



Hi, 



gY J ^\{k)i>\{-k)i> i {-k')^{k'), (3) 



k.k' 



where the bare coupling constant g is negative, to express 
the attractive (s-wave) interaction between the spin j" and 
.j, fermionic atoms. 

After the mean field (MF) approximation and the sub- 
sequent diagonalization of the expression for Hmf, we 
arrive at the following expression for thermodynamic po- 
tential: 



n = 



d 2 k 



d2k (A -4)-^ 



(27T) 
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:i (27T) 



k<ki, k~>k 

k2 d 2 k b 

2 e k> 



d 2 k . , , A 2 



f) 



, (2^) 2 



(4) 



where, for simplicity of notation we have labeled 1= a, 



t= b, and we have defined e 



2 > b k 



Eh 



E,, 



A 2 , and the constant pairing gap is given 



by A = -gf$j7{1>l(-k)1>\(k)) = A*. In Eq. Q 



2 _ k a + ^6 _|_ 1 . , 



±o\/(fc b 2 -^) 2 -16M 2 A 2 (5) 



arc the roots of ejr, and k a = \J2M\i a is the Fermi mo- 
mentum of the species a — a,b. From the equation above 

one sees that the condition for A such that k\^ is real is 

k 2 -b 2 

— AM ■ 

To regulate the ultraviolet divergence associated with 
the first integral in Eq. (|4]) we introduce [3l[ : 



H = H o + H h 



where 



H = J2 4^(k)Mk), 



fc, (7=^,4. 



(1) 



(2) 



is the kinetic (free) part of H and H nt is given below. 
■0j.(fc) and ^a(k) in Eq. ([2]) are the creation and annihi- 
lation operators, respectively, for the f,i particles. To 
assure population imbalance in the system, we have in- 
troduced different chemical potential for the species a 



1 
U 



d 2 k 



1 



(2^) 2 2£ fc + k B r 



(6) 



where U = —g > 0, and es is the 2D two-body bind- 
ing energy. In order to make contact with current ex- 
periments, it is convenient to relate sb to the three 
dimensional scattering length a s . In the scattering of 
atoms confined in the axial direction by a harmonic po- 
tential with characteristic frequency ujl they are related 
by [HIH: 



Cuj, 



-exp 



Z7T 

a s 



(7) 
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w here a s is the 3D s-wavc scattering length, lol = length A parallel to the z-axis [2l| . 
V%r 2 Vb/(A/A 2 ), l L = 1/y/MJE, and C w 0.915. V 

is the amplitude of the periodic potential Vq sin 2 (27rz/A) Thus we obtain an analytical expression for the grand 

generated by two counter-propagating laser beans with canonical thermodynamic potential at T = 0: 



n = n(A,fj,) + e h n(A,h) 

= A 2 



(8) 



In 



/i(V/i 2 + A 2 + M ) -e h 



2/ivV - A 2 - A 2 In 



ft + \//i 2 - A 2 



where n = Q h = 6(h 2 - A 2 ), fi = ^±^, /i = 

tlb ~ fJ - a t and O(x) is the Heaviside step function, defined r i 

as 1 if x > 0, and if x < 0. With these definitions, fi a = ep — . (12) 

Eq. ([5]) may be written as k\ 2 = 2M( l u±v / /i 2 — A 2 ) and 

the condition for real fc 2 2 is now A < h. The value of the free energy at the minimum is 



III. PHASES OF 2D IMBALANCED FERMI 
SYSTEMS AND THEIR STABILITIES 

A. Fixed Chemical Potentials 

In the grand canonical ensemble the chemical poten- 
tials of the two-components are held fixed. This happens 
when the system is connected to reservoirs of species 
a and b such that the particle densities are allowed to 
change. We analyze the ground state of this case for 
balanced and imbalanccd configurations. 



1. Balanced Systems 



n(h = 0,A = A o ) = n o = -n[^ o + L j 1 J , (13) 

where k=M:, whereas the energy of the balanced normal 
state is given by 

0(/i = A = 0) = flfr = —k/Iq. (14) 

A direct comparison between Eqs. ([T3]) and (fT4| shows 
that the supcrfluid state is energetically preferable to the 
normal state for any ^ 0. Since a two-body bound 
state exists even for an arbitrarily small attraction in 
2D [3l| , the pairing instability will always happens in 2D 
balanced two-component Fermi systems at T = 0. 



The gap and number equations are obtained by 
dVl/d A = and n a = —dfl/d/j, a , respectively. For the 
balanced system where h = 0, which implies 6^, = 0, we 
find 



y/}J+A*-H = \e B \, (9) 

and 



v//i 2 + A 2 + A i = 2e F , (10) 

where the two dimensional Fermi energy is defined as 
Ep = ^Sp, with tit = n a + rib, see Appendix A. In the 
balanced configuration n a = rib = n = . Solving these 
two equations self-consistently we arrive at 



A = y/2e F \e B \, (11) 

and 



2. BCS-BEC Crossover in Balanced Systems 

The two well known [3l], HH equations (fTT|) and (fT2|) 

reveal that the BCS-BEC crossover can be accessed by 
varying cb from weak to strong interaction regimes. The 
BCS state is reached in the weak attraction limit or high 
density (ep), \cb\ << tp, where no = ep and the energy 
gap is found to be 

A = 2^e F E A e~^j << £Fi (15) 

where Ea = , and A is a momentum cutoff. To obtain 
Eq. (fT5l) from Eq. (fTTI) we have made use of Eq. ©. The 
opposite limit, of very strong attraction or low density, 
| 6b I >> Ep, results in the formation of composite bosons 
with fiQ = — | eg |/2. Since the BCS and BEC regions are 
characterized by positive and negative chemical poten- 
tials, respectively, a clear distinction between these two 
regimes is the value of the binding energy at which /i 
changes sign i.e., |es| = 2ep. 
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3. Imbalanced Systems 

We now turn our attention to the cases where h ^ 0. 
The free energy of the imbalanced normal state. 0(/i, A = 
0) = Of, is found to be 

In the limit A^Owc find k x (A = 0) = y/2M {p - h) = 
k a F < k F and fc 2 (A = 0) = y/2M(fi + h) = k b F > k F , 
where kp = \J2M\i. We are considering fi a and h posi- 
tive. Thus the equation above is written as: 

°?-L <^ + i (17) 

which gives, after the integration in k, the free energy 
of a (normal) two species gas of fermionic atoms in two 
dimensions 

nf^.w) = -|[( M -M 2 + (/i + ^) 2 ] 

= -f[M 2 +M 2 ]. (18) 
The gap and number equation now read 

in ( ^J^-A + Qh]n ( ^^ZF] = o, 
\ \e B \ J \h-Vh 2 -A 2 J 

(19) 



Vm 2 +A 2 +m = 2£f. (20) 

Note that for h < A Eq. JTS} is reduced to Eq. ©. 
Seeking now solutions where h > A we then have to 
solve 

V/x 2 + A 2 - /i = ft - Vfe 2 - A 2 

|ea| h + Vh 2 - A 2 ' 1 J 

Solving these equations self-consistently we find 

A s (h) = ^A (2h- A ). (22) 

From the equation above we see that the "Sarma state" 
(as explained below) gap and the chemical potential im- 
balance exist in the ranges 

< A s (h) < A , (23) 



~Y<h<A , (24) 



where the upper limit is imposed by the existence of real 
k 2 

From the graphical analysis of f2 as a function of A 
for various asymmetries, shown in Fig. [TJ one sees that 
the first curve, from top to bottom, is the usual balanced 
system with its minimum at Ao. Increasing the imbal- 
ance h and keeping /z positive (to guarantee that we are 
in the BCS regime), the minimum is still located at Ao 
up to a maximum or critical imbalance h c , after which 
there is a quantum phase transition to the normal state 
with A = 0. h c is find through the equality fio = Of, 
which yields: 



hl = f ^\e B \ + (^y , (25) 

from which one easily finds plugging in the equation 
above fj, from Eq. (IT21 in the BCS limit (|e B | << e^): 
h - 

We have seen above that, as happens in 3D [H, Q, 
the BCS phase turns from stable, while h < h c , to 
metastable, when h > h c , in which case the normal phase 
is stable. Besides these two phases, there is an unstable 
phase, known as Sarma state, corresponding to a local 
maximum of fl versus A, that is located between the 
BCS minimum at Aq and the normal phase with A = 0. 



-1.9 



Q 




-2,3-1 , , , , , 1 , 1 

0,2 0,4 0,6 0,8 

A 

FIG. 1: The thermodynamic potential $1 as a function of A 
for various chemical potential asymmetries. The top curve 
is for h — 0. The imbalance is being increased from top 
to bottom and in all curves for h < h c the minimum is at 
Ao. The curve for h = h c has two minima. After h c , at 
which there is a quantum phase transition for the normal 
state, f2(A = 0, h > h ) < fi(A , h < h c ). 
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B. Fixed Number Densities 

In the previous subsection we discussed the stability 
of several phases with fixed chemical potentials [i a and 
fib- We investigate now the situations where the num- 
ber densities n a and rib are fixed, since it is appropriate 
to describe trapped gases in current experiments. It is 
convenient to introduce now the "magnetization" defined 
(see Appendix A) as m = — = n& — n a , which will 
be needed to describe the energy of the Sarma phase in 
terms of the number densities, as discussed below: 



m= — Vh 2 -A 2 . (26) 

Solving equations ([20)) , (|2"Tj) and (|26|) sclf-consistently, we 
find the gap As(m) and the "average" chemical potential 
fi s of the Sarma phase 

A s (m) = ^Al-^p-m, (27) 

where Ao, the gap parameter of the balanced system, is 
given by Eq. (fTTj) . and 



(28) 
I 



Note that the equations above reduce to Eqs. (fTTj) and 
([T2"|) . respectively, since in the limit m — > 0, As — > A , 
as it should be. From Eq. ([27]) it is easy to obtain the 
windows for the Sarma gap and m: 



< A s (m) < A , (29) 



< m < m max , (30) 

where m max = is the maximum value for the den- 

sity imbalance in the Sarma phase, which is easily ob- 
tained from Eq. ([27]). 

The expression for the energy E (not the free energy 
or thermodynamic potential) of the homogeneous Sarma 
phase is obtained as: 



E(n T , m, A) = fl(fx, h, A) + fi a n a + fi b n b . (31) 

To write E = E{tit, rn), we express fi a n a + = /i«T + 
hm, with jx = ii{nr), h = h(m) to find 



E(riT,m,A) = E s (n a , rib) — Q(tit, m, A) + \itit + hm 



™ A 2 



Ae F J 



A 2 

-2e F [e F --^ 
4e F 



- 2 



A| m + A| In 



yjfh 2 + 



(32) 



where m = ^jj- . 

To investigate which one is the ground state of an 
imbalanced system with different densities n a and rib, 
we have to see which energy of the following states we 
are considering is smaller, the normal E N (n a ,rib), see 
Eq. (|36|) below, the homogeneous Sarma E s (n a ,rib), or 
the phase separation (PS) E PS {n a ,rib) state. The PS is 
an inhomogeneous phase where given n a and rib densi- 
ties in a trap, a fraction 1 — x of the 2D (real) space is in 
the BCS phase where both species have a common den- 
sity = n^ cs = n and the rest of particles are in 
the normal phase occupying the fraction x around 1 the 



1 A surface energy is related to the interface between the normal 



BCS core, with densities n a and fib- The most favored 
composition minimizes: 

E PS = Min xfi [(l - x)E(n) + xE N (fi a , n b )\, (33) 

where the space fraction is obviously limited to the range 
< x < 1. The energy of the balanced superfluid phase 
is found as 



and superfluid phases [3414371 ] . The surface tension associated 
with this interface has been interpreted as responsible for defor- 
mations in highly elongated small samples [14t llal . These results 
are not consistent with the observations reported in Ref. [Itl . We 
neglect the surface energy in the present analysis. We intent to 
come to this issue again soon taking into account corrections 
beyond the mean-field approximation. 
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Phase Competition 



E(n) = fi(A, h = 0, n) + \in T 



^{a 2 
MA 2 



4tt 







Kl) 


- 1 





2ei 



n(2n) 2 
2M 



Since E(n) which enters E PS has to be written at its 
minimum, the equation above turns out to be: 



There are some limiting cases where the comparison 
between the phase separation and Sarma states can be 
done analytically and exactly. The first one is obtained 
when m = 0, giving E s (n) = E PS {n) = E BCS (n). The 
(34) second case happens when to = m max , which reduces 
the Sarma phase to the normal phase with A = 0. This 
is easily seen from Eqs. (|30|l and (|2l?j) . Then we define 
AE = E PS - E N , which is written as 



E(n) 



TTTlT 

2M 



47T 



\ZB\n, 



(35) 



where the gap parameter of the BCS phase is given by 
Eq. (fTTj) and written as A (n) = 2^/-jj|es|n. 

The energy of the normal state entering in Eq. ([55)1 is 
easy to find from Eq. (|18|) as 



E (n a ,n b ) 



Note that since 



fl N (n a ,n b ) 
1 

2k 



K + ni). 



(36) 



and n b = 
pJV _ t™? 



, the equa- 
tion above for to = gives E N = agreeing with the 
result of Eq. ((MJ) in the limit A ->• 0. 

As we mentioned before, the n a and n b particle densi- 
ties are accommodated in the trap as 



n a 
rib 



Xfi a + (1 — x)fl, 

xn b + (1 — x)h. 



(37) 



Then we rewrite Eq. (|36p as 



E (n a ,h b ) 



1 

2^ 



n a - (1 - x)n 



Tib — (1 — x)n 



(38) 



It is interesting to note that Eqs. (|3"5"T) and (|3l)|) arc 
easily obtained from Eq. (|32p in the limits m — > and 
to — > m max , respectively. 

In the minimization of E PS we shall find i « and 
n ~ n a ( or nt) if n b n a , which means that the entire 
system is in the supcrfluid phase. In the other limit, if 
n b » n a , the system will be in the normal state, so we 
shall find in the minimization of E PS x « 1 and n w 0. 
In the intermediate cases the system will phase separate, 
i.e. if the imbalance is less than the critical value, as 
experiments have shown for 3D gases. 



AE 



— - E(n) 

K 



(x-1) 



Minimization with respect to x yields: 



/ (tit - n) 2 - 2n a n b + h 2 . . 

Xmm "V 2(n 2 - KE(h)) ' 1 ' 



An upper bound on AE can be obtained by setting 
the density of the BCS component of the PS state n = 
n a and writing n b = n a + m, where for low imbalance, 
to < n a ,n b , Eq. ([ID)) gives: 



and 



AE(x m i n ^ 



IP 



^ MA 2 (n a ) 



(41) 



3 „2 



4tt 



3tt 
2M 1 



(1 



,) 2 < 0. 



(42) 



Unlike the three-dimensional case where the result ob- 
tained is an approximation, valid only for small asym- 
metries to << n a ,rih, where terms of order TO 3 /na 4 / 3 
have been neglected d, Q , the result above is exact and 
valid for any to < n ai n b . These results agree with the 
numerical analysis made to construct the phase diagram 
depicted in Fig. ©. 

Let us now show the results of the minimization of 
Eq (|33[) with respect to x and n for a given n a and n b . 
The results are shown in Fig. ([2]). 

We show in Fig. (J3j) an illustration of what may 
be found experimentally in an imbalanced 2D gas of 
femionic atoms. The first picture represents the balanced 
system with to = 0. The following pictures represent a 
set of experiments with fixed es/e^ = 0.0318 with in- 
creasing imbalance. From the second picture to the fifth 
the system phase separates and the sixth picture is for 
the critical imbalance m c representing the normal phase. 



FIG. 2: The phase diagram of a zero temperature imbalanced 
gas of fermionic atoms in 2D. The vertical axis |es|/e_p with 
m/riT = 0, is the (balanced) BCS phase, whereas the hori- 
zontal line |es|/eF = represents the (imbalanced) normal 
phase. Note that for each value of |es|/eF phase separation is 
the ground state until m c . Above this value phase separation 
is not energetically favorable, and the normal phase has lower 
energy. 



2. BCS-BEC Crossover in Imbalanced Systems 

Now we investigate the BCS-BEC crossover in an im- 
balanced gas of fermionic atoms. As we have seen, there 
are three possible stable phases with fixed number densi- 
ties. For zero imbalance between the densities, the only 
phase present is BCS (for any value of the interaction 
between the spin up and down fermionic atoms) and the 
crossover is that discussed in Subsection MI A 21 which is 
governed only by cb- For imbalances where phase sepa- 
ration is possible, the BCS-BEC crossover can be realized 
in the superfluid core of the separated phases, again vary- 
ing 6b in fiQ = ^jj- — , where h is the BCS density that 
minimizes the PS formation, as described above Eq. (|33|) . 
For large enough asymmetries, phase separation is not 
supported and the entire system is in the normal phase 
with no possibility of pairing formation and the crossover 
to the BEC regime. Then, when m ^ the BCS-BEC 
crossover can not be accessed only by varying the biding 
energy since the imbalance plays an important role now. 



IV. CONCLUSIONS 

We have investigated superfluidity in two-dimensional 
imbalanced Fermi gases. With exact expressions, we have 
constructed the zero temperature mean-field phase di- 
agram of these systems. We show that in the grand 
canonical ensemble, BCS is the ground state up to a 



m/n^O.08 m/n^O.12 




m/n^O.16 m/r> T = 0. 179 



FIG. 3: Evolution of the phases of an imbalanced gas with 
fixed |es|/eF = 0.0318 as a function of increasing m/nr- 
The first picture, with m/nr = represent the BCS phase. 
Increasing m/nr the system phase separates (pictures 2 to 5), 
and increasing the imbalance even more, the system finally 
goes to the normal phase at the critical imbalance m c /riT = 
0.179, as seen in the last picture. 



critical chemical potential imbalance h c at which there 
is a first-order phase transition to the normal state. In 
the canonical ensemble, relevant to current experiments, 
we also show that, as found theoretically and experimen- 
tally in 3D, for a fixed two-body binding energy es, as 
the imbalance m is increased from zero the stable states 
are BCS, phase separation and normal phase. For bal- 
anced systems, the BCS-BEC crossover is governed only 
by es- However, for imbalanced systems, the BCS-BEC 
crossover will depend on es and the imbalance m. For a 
fix value of e^, the crossover is possible only for asym- 
metries supported by phase separation. Above a critical 
value of the number densities imbalance m Cl phase sepa- 
ration is destroyed in favor of the normal phase, as well 
as the path to BEC regime. Thus, for a given m < m c 
such that the system phase separates into BCS and nor- 
mal phases, the BEC regime can be accessed by vary- 
ing es from the weak to strong couplings. A smaller 
and smaller circular core shall be visibly apparent, as es 
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is increased, signalizing the condensation of the Cooper 
pairs. This is a direct prediction of our results that could 
be tested experimentally. Another interesting phenom- 
ena that could also be verified experimentally in a 2D 
trapped imbalanced Fermi gas [H| , emerges in the oppo- 
site range, that of very large imbalance m > m c or, in the 
unities of Fig. [21 m/nx ~ 1. In this situation, the minor- 
ity species a correspond to impurities imcrscd in the ma- 
jority species b Fermi liquid. As the interaction |es|/e_F 
is increased, the impurities a surround themselves with 
a localized cloud of environment b atoms, forming Fermi 
polarons (39| . 



an, on , on , A 

o/ia o^i b OA 

Since the densities have to evaluated at the minimum 
of Q, i.e., at A , the last term vanishes. The particle 
densities are given by n a = — Jjp-. Then, using \x a = 
H — h and fii, = /i + h in Eq. (|43|) we find 
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VI. APPENDIX A: DERIVATION OF THE 
MAGNETIZATION AND TOTAL NUMBER 
EQUATIONS 



n T = n b + n a 



d/i 



(45) 



The total derivative of = f2(/i a , A*6, A) is written as 
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